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Abstract 



The four observables associated with gravitational lensing of distant quasars by in- 

<N . 

O . tervening galaxies: image splittings, relative amplifications, time delays, and optical 

°. 

depths, provide separate measures of the strength of the gravitational constant G 
\ at cosmological distances. These allow one, in principle, to factor out unknown 



lensing parameters to directly to probe the variation of G over cosmological time. 



We estimate constraints on G which may be derivable by this method both now 
and in the future. The limits one may obtain can compete or exceed other direct 



limits on G today, but unfortunately extracting this information, is not independent 



of the effort to fix other cosmological parameters such as H and Qq from lensing 
observations. 
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1. Introduction 

The gravitational constant, G, is the poorest measured fundamental con- 
stant in nature. In fact, it may not even be a constant at all. The exceedingly 
small value of G, coupled with the large value of the age of the universe encour- 
aged speculation early on, first following Dirac and then spurred by the advent 
of Brans-Dicke cosmology, that the two quantities may be somehow tied together 



[ Pirac 1937| , Pyson|| . Moreover, because classical general relativity cannot be quan- 
tized, there has been a recurring interest in the possibility that GR arises as the 
low energy limit of a more fundamental theory. In such a theory, the gravitational 
constant may arise dynamically, associated with the vacuum expectation value of 
some field (or dynamics of some internal space). Since this dynamical value may be 
time dependent, so may G. Over the past year, largely as a result of considerations 
based on extensions of the original old inflationary models ||Guth 198T], |Linde 1982 



Albrecht and Steinhardt 1982| , |Steinhardt and Accetta 1990|| , there has been a re- 



newed interest | |La and Steinhardt 1989| | in the possibility that the gravitational 



constant has varied on cosmological timescales. 

There exist several sensitive direct probes of a monotonic change in the grav- 
itational constant during the present epoch, including the use of pulsar timing mea- 
surements and radar experiments, all of which suggest that G/GH < 0.4 today. 
RShapiro 1964] , [Shapiro et al. 1971| , Polling 1987] , |R easen 19831 , Pamour et al. 1988| 



At the opposite extreme, calculations of primordial nucleosynthesis put indirect 
limits on G during the first seconds of the big bang expansion from limits on the 
observed Helium abundance jAccetta et al. 1990|] . If the variation of G has followed 



a constant power law in time, the latter limit (G/GH < 0.01) is stronger than the 



direct limits on the variation today . What has been lacking however is any way to 
directly probe the value of G at times between these two epochs. Since it has even 
been proposed that G may oscillate in time [Accetta and Steinhardt 1991| , a direct 



measure of G at intermediate times would be of great interest. It is the purpose 
of this paper to suggest that observations of gravitational lensing could, in princi- 
ple, provide such a measure, and to investigate the realistic limits which it may be 
possible to obtain. 

On first thought it is not clear that lensing can constrain G. While the 
bend angle which light rays are subject to is directly related to the strength of the 
gravitational constant at the time light rays pass the lensing object, the quantity 
which enters into all formulas is the product GM, where M is the mass of the 
lensing object. Unless M can be determined independently a separate extraction 
of G seems impossible. However, it is not the actual bend angle which is directly 
observed in gravitational lensing. All lensing observables depend also (in somewhat 
different ways) upon the distance of the lensing galaxy and the quasar as inferred 
from their redshifts. The distance redshift relation depends upon the time-averaged 
value of G, which for redshifts of 0(1) can be a significant fraction of the lifetime 
of the universe. Thus for any lensing system a prediction of one observable based 
on a measurement of another can give a signal of the time variability of G. What 
remains to be seen however is exactly how sensitive such a comparison is, and how 
much it depends on our knowledge, or lack thereof, of cosmological parameters such 
as the Hubble constant H , the density parameter Q Q , and even the cosmological 
constant A. 

The organization of the paper is as follows: in section 2 we outline our nota- 



tions and conventions and introduce the models we will use. In section 3 we consider 
constraints from lensing statistics and in section 4 we discuss constraints which can 
be derived from individual lensing systems. Section 5 contains our conclusions. 

2. Cosmology and Lens Models 

The observables of interest in gravitational lensing depend upon the combi- 
nation GM (where M is the mass of the lensing galaxy) and the distance to the 
galaxy and source. If it is assumed that the bending occurs predominantly as the 
light rays pass through the local region of the lensing galaxy and thus the time 
required is much shorter than the time scale over which G varies significantly, the 
effect of the variation of G will be to replace GM by GiM (where Gi is G at the time 
of lensing) and also to alter the distance-red shift relation. Measuring the first effect 
is cleanest, in principle, because it is not dependent upon cosmo logical modelling. 
Unfortunately, unless there is an independent way to determine the mass of the 
lensing galaxy the first effect alone is unmeasurable. Since both velocity dispersion, 
and stellar luminosity will also depend upon G, there are no observables which seem 
to allow M to be independently extracted. 

Hence, to proceed, we must consider some specific cosmological model, incor- 
porating a variable G. We will consider for definiteness a Brans-Dicke (BD) theory 
(this is perhaps the simplest viable extension of GR with a varying gravitational con- 
stant and is often used in connection with extended inflationary models). While our 
discussion will be in terms of BD cosmology, the general features should be character- 
istic of any model with varying G. In particular, these ideas could be applied to any 
theory based on the Friedman-Robertson- Walker metric (which has gained more ex- 



4 



perimental support recently from the isotropy of the cosmic microwave background) 
with an evolution equation for the scale factor determined by the equation of state 
of matter which also incorporates a varying value of G consistently in the equations 



of motion (e.g. |pyson| , [Beckenstein 1977[ [Beckenstein and Meisels 1980 



In the BD cosmology the line element is the usual FRW metric 

ds 2 = -dt 2 + R 2 {t)[d X 2 + Skixfdfl] (1) 

where Sk = sinh(x), x, sin(x) for k = — 1, 0, 1. Einstein's equations are modified and 
a new dynamical field 0, with G ~ <p~ l , is introduced. For large time the general so- 
lution of the Brans-Dicke field equations will be matter dominated and in many cases 
of interest (e.g. ||Weinberg 1972| , [La and Steinhardt 1989| , |Steinhardt and Accetta 1990|| ) 



<p ~ R a , where a will be a function of the Brans-Dicke scalar-tensor coupling 
constant 3 , ou, which tends to zero in the limit uj —>■ oo (where Einstein's theory 
is recovered). 

| 3 The independent limits on the scalar-tensor coupling constant in the simplest 
Brans-Dicke theory are already far more stringent than we will place from the vari- 
ation of G [[Reasonberg 19791 , but our purpose here is to use this model merely as 



a testing ground to explore the sensitivity of lensing parameters to G.\ 

In the (k = 0) examples cited above a = (1 + u)~ l . The evolution equation for the 

scale factor in a matter dominated epoch is 

(R\ 2 k 8rrG f2u + 3\ fR \ 3+(7 ,a 2 u . ( R, 



If we define r\ = R/Rq and 



we can rewrite this as 



V 2 + (O - l)Hl = Q H 2 oV - {1+a) (4) 



2 N 



^ = ^o-l)H^l + a-^j (5) 

Note that if one instead were simply to allow G to vary as a power law and 
were to use Einstein's equations unchanged one would obtain the same result (f|.[5|). 
except without the last factor on the r.h.s. of ([5]). For a fixed value of Qq and Ho, 
it is this factor which causes Ro to vary with ui. 

The measure of distance we will use is the angular diameter distance 

dA = y^-s k (x) (6) 

which assumes that the lensed rays traverse a mean filled "beam" . (Similar, but alge- 
braically more complex constraints can be obtained in the case of an "empty" beam 
approximation, where affine angles and distances are used (i.e. see ||l'urncr et al. 198~i . 



[Krauss and White 1991|1 ). The mean filled beam approximation is probably closer 



to the actual situation, however, it has been shown that the uncertainty due to 
dumpiness of matter in the beam trajectory can be one of the main sources of un- 



certainty in the analysis of individual lensing systems fAlcock and Anderson 1985 
If the Hubble constant were independently measured it is possible that this uncer- 
tainty could be reduced since it also enters into the determinations of H from time 
delays in lensing systems. 

Once the evolution equation for the scale factor, R(t), is specified we can 
solve for the distance redshift relation in the usual way [[Weinberg 1972|| . For the 



special case k = we obtain a simple expression for the distance as a function of 
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x = 1 + z, (similar expressions for Oo ^ 1 can also be obtained.) 

where f3 — (1 + u)/2. The distance is plotted in figure 1 for (3 — 0.4,0.5,0.6, 
and reduces to the usual expression ||Turner et al. 1984| in the limit of constant G 



{uj — > oo, (3 — > 1/2). As a guide to the expected magnitude of (3 one would like to 
obtain sensitivity to, notice that an assumed variation-since-lensing of 

^ = ^ = 20 %^ =**H~0.2 (8, 
G G log xi 

for lenses z\ ~ 1.5. This corresponds to /3 ~ 0.6. If we take the age of the universe 
to be t = 2/(3 + (j)H l = 10 10 yr, this then gives 

G/G\ = -aH ~ 10- n yr _1 (9) 

which is comparable with other direct measures of G/G. 

We will consider two simplified lens models in what follows: the point mass 
and the isothermal sphere lenses. The point mass lens is chosen for its simplicity, the 
isothermal lens because the flatness of rotation curves of galaxies suggest p ~ r~ 2 is 
a reasonable approximation to galactic mass distributions (at least asymptotically). 
While for any actual lens system these models are overly simplistic they serve to 
illustrate the main points. The observables for these lens systems which we would 
want to examine for sensitivity to (3 are: time delays between images, angles between 
images and ratio in brightness of the images, as well as the (differential and total) 
optical depth for lensing. 
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3. Lensing statistics 

The formalism appropriate to gravitational lensing statistics was first devel- 
oped in fTurner et al. 1984]] and latter generalized to arbitrary Robertson- Walker 



cosmologies in | |Gott et al. 1989| (see also [[Krauss and White 1991| | for a recent pre- 



sentation). The key quantity is the optical depth, or integrated probability of 
lensing, r, assuming a non-evolving population of galaxies, modelled as singular 
isothermal spheres. This depth is relatively free of matter clustering uncertainties 
Alcock and Anderson 198~5f but is sensitive to variations in the distance-redshift 



relation, which makes it a good probe of cosmology. 

As an example consider the expression for r, for a k = universe. Including 
the ft dependence of the bend angle (a ~ G) we obtain 



r (y = 1 + z) = — \ y dx x'^^x 13 - l) 2 



ft 2 J i 



y 13 — x 13 
~yT^T 



2 



F = 71o7TQ!q-Rq (10) 



where H R = 1 and n is the comoving number density of galactic lenses, which 
we assume in this instance are all identical (non-evolving) isothermal spheres pro- 
ducing identical bend angles oiq. The bend angle «o is related to measured velocity 
dispersions of nearby galaxies today. For further details of these definitions see 
[ Turner et al. 1984 1 . Equation (|lO"l) reduces to eqn (2.26c) of [ Turner et al. 1984 



m 



the limit /3 — > 1/2. The integral can be done analytically but the result is cumber- 
some and is not shown here. The optical depth vs redshift is shown in figure 2 for 
ft = 0.4, 0.5, 0.6. As can be seen the variation with ft is slight making this a poor 
measure of G/G. A similarly small dependence on ft is shown by the differential 
optical depth dr/dz. Since the distance-redshift relation becomes less ft dependent 
as flo decreases we expect the variation in r to be less than above when Qq < 1, al- 
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though this is somewhat offset by the (3 dependence of F coming through Rq. Thus 
variations of G going as a power law in time have little effect on lensing statistics, at 
least at the level where these statistics are likely to be determined in the forseeable 
future. 

4. Individual systems 

A better hope of constraining G comes from examining the observables asso- 
ciated with multiply imaged quasars (i.e. see [Hewitt et al. 1988|| ). Specifically we 



will be interested in the observables: time delay, image splitting and image magnifi- 
cation for our two model lenses. The strategy will be the following: each observable 
will depend both on GM, and on <1a{G). If we have more than two observables 
for each system, then we hope to overly constrain the system so that we can check 
for consistency between the different determinations of these quantities from each 
observable. 

1) Point mass lens. 

For these lenses the time delay At and ratio in magnitude of images r are 
related through [[Krauss and Small 1991| | 



At = 2(73/(1 -i ;)[(/■-— )+k>g(r)] (11) 



(the (1 + z) factor is absent in microlensing fKrauss and Small 1991|| ) so these two 



parameters can be used to infer GM at the time of lensing, independently of the 
cosmological distance-redshift relation. If it is not possible to measure At in the lens 
system, or if the measure has a large uncertainty, GM must be obtained some other 
way, e.g. from virial velocity measurements. Any limit on f3 will depend on how 
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well this quantity is known. 

Given GM we can use the observed angular splitting of images and the 
relation 

. „ 4GM r - 1 

c 2 S ,^ r i/2( r _ 1 )_ 2r 

to determine S = Ds/Dls where Ds and Dls are the angular diameter distances 

from the observer to the source and from the lens to the source respectively. This 

is a function only of the (known) redshifts, Q and /3, e.g. in the k = case 

S=T^L = 4^4 (13) 



so(xs - Xl) 



X L X S 



so a knowledge of the redshifts allows a determination of f3 (up to clumping un- 
certainties l|Alcock and Anderson 1985)1 ) if we assume a value for Qq (or conversely 



a determination of Qq if we know (3). Notice that S is a ratio of distances and so 
is independent of H . As an example if we take Zl — 1, z$ — 3 then for k = 0, 
mean filled beam, S is a monotonically increasing function of (3 varying from 2.3 to 
2.5 as (3 runs from 0.4 to 0.6 as can be seen in figure 3. Given the above and the 
fact that typical image splittings can be ~ 3" — 7" it is not impossible that a good 
measurement of the angular splitting could limit (3 to be in the range competitive 
with other direct probes of G. 

2) Isothermal Sphere 

For the somewhat more realistic, isothermal sphere model the situation is 
simpler (in principle). If the velocity dispersion of the lensing galaxy, mi, is known, 
say from measurements of the rotation curves, a measure of the angular splitting 
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allows us to immediately infer S: 



2a 87Rrjj „ 8vr /a N \ 2 



A9 = ^ = ^ s = MfJ < 14 > 

In fact a simultaneous measurement of the time delay, which gives us D = D^/S 
(where Dl is the angular diameter distance from the observer to the lens), 

At = 32tt 2 (1 + z) f^-) — (15) 



and the image splitting would in principle allow us to measure both (3 and VLq (up 
to uncertainties in H , because Dl has the dimensions of distance and hence is de- 
pendent upon H ) because the dependence of D and 5* on (3 and Qq is different 4 . 
We expect the strongest constraint on (3 for fixed Qq to come from S however. 
| 4 See figure 4.| 



The singular isothermal sphere (SIS) model is probably still too naive to 
apply to actual individual lens systems. One should at least include the effects of 
a finite galactic core [[Hinshaw and Krauss 19871 , IKrauss and White 1991| . Alterna- 



tively, a more complicated, but more general model, the 'elliptical lens' 
[Narayan and Grossman 1988|| , is available for use in extracting these quantities for 



individual galactic lenses. Nevertheless, the SIS model should give a general idea 
of the methodology to be used, and the possible sensitivity to (3. Using these other 
models in an application of these ideas to actual lenses would merely require replac- 
ing the above equations for S and At with somewhat more complicated equations 
which would include the lensing parameters fit by the observations. 

Because of the simplicity of the SIS model we only had to make due with 2 
lensing observables to overconstrain the system. We note however that the ratio of 
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image amplifications itself is also dependent on S and (3, and so can also be used to 
probe for consistency when more complicated fits to galactic lenses are required. 

5. Conclusions 

While we have demonstrated here that gravitational lensing provides in prin- 
ciple a direct sensitivity to variations in G over cosmological time, our results suggest 
that to be competitive with limits on G at the present time, lensing parameters must 
be extracted from observations at the level of 10% or better-a daunting but not im- 
possible task. Statistical measures such as the optical depth do not seem sufficiently 
sensitive to G/G, the effects of reasonable changes in G being swamped by larger 
uncertainties from our present lack of knowledge of Qo and H . For individual lens 
systems a simplified model suggests that it may be possible to see variations of the 
order G/G ~ lCr n /yr if accurate measures of the angular splitting, amplifications, 
and perhaps also time delays become available for a system with source and lens at 
relatively high redshifts (z s ~ 2, 3 and zl ~ 1)- Pessimistically, it is worth noting 
that a possible variation in G is yet one more uncertainty which could limit one's 
ability to extract H Q from measurements of time delays in individual systems. On 
the other hand, if H and Q Q are measured reliably by independent means, one's 
ability to probe for variations in G will improve. 

Nevertheless, in spite of the limitations of this method, it is worth empha- 
sizing that it does provide perhaps the only 'direct' probe of variations in G during 
intermediate times between the present epoch, and the nucleosynthesis era in the 
very early universe. We have placed 'direct' in quotation marks because as we have 
demonstrated, one's ability to extract information on G is intertwined with our 
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knowledge (or ignorance) of the proper cosmological model for the evolution of the 
universe during this time. In this regard, we also note that for our analysis, we 
used as an analytic tool to probe the sensitivity of lensing, a simple Brans-Dicke 
cosmological model, which in fact is already ruled out by other constraints for the 
parameter range which would produce the level of time variation probed here. In 
this case, G would vary as a simple power law with time. We expect our results 
would be applicable for any similar, more viable, model. Of course, there are other 
possibilities, including an oscillatory behaviour of G with a cosmologically interest- 
ing period (i.e. [Accetta and Steinhardt 1991|| ). One would need specific models to 
perform an analysis similar to that performed here, but it may be that for such 
scenarios, gravitational lensing could provide sensitive limits, by comparing results 
obtained from lensing systems at different redshifts. 
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